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! Abstract 

H ■ 

We present a method for estimating the edge of a two-dimensional bounded set, given a fi- 
nite random set of points drawn from the interior. The estimator is based both on projections 
on C 1 bases and on extreme points of the point process. We give conditions on the Dirichlct's 
kernel associated to the C 1 bases for various kinds of convergence and asymptotic normality. 
We propose a method for reducing the negative bias and illustrate it by a simulation. 
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1 Introduction 



We address the problem of estimating a bounded set S of R 2 given a finite random set £ of 
points drawn from the interior. This kind of problem arises in various frameworks such as classi- 
fication [10], image processing [14] or econometrics problems [3]. A lot of different solutions were 
ly-) | proposed since [6] and [16] depending on the properties of the observed random set £ and of the 
unknown set S. In this paper, we focus on the special case where E is the set of points of an 
homogeneous Poisson process whose support is S = {(x, y) G R 2 | < x < 1 ; < y < f(x)}, 
where / is an unknown function. Thus, the estimation of subset S reduces to the estimation of 
function /. Let us note that this kind of support was already considered in [6]. In the wide range 
of nonparametric functional estimators [2], piecewise polynomials have been especially studied 
^ [13, 14, 15, 5, 9] and their asymptotic optimality was established under different regularity con- 

ditions on /. 

The first support estimator based on orthogonal series appears in [1]. Its properties are exten- 
sively studied in [12] in the case of Haar and C 1 bases. The expansion coefficients estimation 
requires the knowledge of the process intensity. This is a limitation which is avoided in [7] in the 
case of the Haar basis by considering a coefficient estimation based on the extreme points of the 
sample. In this paper, a similar study is carried out in the case of C 1 bases. The estimator can 
be written as a linear combination of extreme values involving the Dirichlet's kernel of the C 1 
basis. A close study of the extreme values stochastic properties as well as a precise control of 
the Dirichlet's kernel behavior allow one to establish general conditions for various convergences 
and asymptotic normality of the estimator. Our results are illustrated for the trigonometric basis 
case. Note that the model proposed here is well-adapted to the estimation of a bounded star- 
shaped subset of the plane. Let D be such a domain. Then, there exists a convex subset C of 
this domain, called the kernel of D, from which the whole boundary dD of D can be seen. If 
we assume an interior point of C to be known, the use of polar coordinates allows to reduce the 
problem of estimating dD to the problem considered here, that is the estimation of a function 
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/, with the particularity that /(O) = /(l). In such a situation, extreme points observed in the 
neighborhood of x = bring information on the behavior of / in the neighborhood of x = 1 and 
vice versa. Thus, an estimation relying on the trigonometric basis is specially well adapted. 
This paper is organized as follows. Section 2 is devoted to the definition of the estimator and Sec- 
tion 3 presents some basic results on extreme values and Dirichlet's kernels. The mean integrated 
square convergence of the estimate is briefly studied in Section 4 and the asymptotic normality is 
established in Section 5. In Section 6 a very simple bias correction is proposed, illustrated in [8] 
by a simulation. 



2 Definition of the estimator 
2.1 Preliminaries 

Let N be a Poisson process with a mean measure fx = cA, where the intensity parameter c is 
unknown, A is the Lebesgue measure, and the support of N is given by: 

S = {(x,y) eR 2 | 0<x< 1 ; 0<y</(x)}. (1) 

We assume that / is measurable and satisfies 

< m = inf / < M = sup / < +oo, (2) 

[0,1] [0,1] 

which entails the square integrability of / on [0, 1]. In the sequel, we will introduce extra hypoth- 
esis on / as needed. 

Let (ej)j g N be an orthonormal basis of L 2 ([0, 1]). The expansion of / with respect to the basis is 
supposed to be both L 2 and pointwise convergent to / on [0, 1]: 

+oo 

Vx€ [0,1], f{x) = ]Ta iej (x), (3) 

i=0 

with 

Vz > 0, oi = / ei{t)f{t)dt. (4) 
J o 

We denote by K n the Dirichlet's kernel associated to the orthonormal basis (ej)jgN defined by 

K n (x,y) = '$2e i (x)e i (y), (x, y) £ [0, l] 2 , (5) 

i=0 

where (h n ) is an increasing sequence of integers. The trigonometric basis will provide us with an 
important example to illustrate our convergence results. It is defined by 

eo(x) = 1, e2fc_i(x) = \/2cos2A;7rx, e2fe(x) = v / 2sin2A;7rx, k > 1, (6) 

and we shall suppose for convenience that h n is even. This leads to 

sin(l + h n )ir(x - y) 



K n (x,y) 



sin-7r(x — y) ^' (7) 

1 + K x = y. 



The speed at which the sequence (a?;) decreases to is linked to the regularity of /. In the case 
of the trigonometric basis, if / is a function of class C 2 then 



O (k- 2 ) , (8) 



(see [4]). The estimator is built in two steps. First, in subsection 2.2, / is approximated by 
a sequence (/ n ) obtained from its expansion with respect to the orthogonal basis. Then, in 
subsection 2.3, an estimator f n of f n is proposed. 
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2.2 Approximation of / 

Let (k n ) be an increasing sequence of non-negative integers such that k n = o{n). Divide S into 
k„ cells D„ r where: 



D n ,r — { i.X,y) G S | X G In,r }j Ai,r 



r — 1 r 



Each coefficient a« is approximated by discretizing (4) according to: 

h. 

^ pJt.A\(D^ A rr„ = - 

2fc 



, r=l,...,/c n . (9) 



2r — 1 

a i,fcn = e i( x r)K D n,r), X r = — . (10) 



Then, the expansion (3) is truncated to the h n first terms leading to 

fn{x) = a itkn ei{x), x G [0,1], (11) 

i<h n 

which can be written in terms of the Dirichlet's kernel: 

f n {x) = J2 Kn(Xr, x)X(D n>r ), X G [0, 1]. (12) 

r=l 

Let us emphasize that the approximation f n of / only depends on the basis (ej)j G N through its 
Dirichlet's kernel. The next step towards the definition of the estimator consists of estimating 

A(A», r ). 

2.3 Estimation of f n 

Let N* n denote the superposition N\ H + N n of n independent copies of the point process N, 

and Ti n the set of points generated by N* n . For r = 1, . . . , k n , consider the maximum X* r of the 
second coordinates of the set of points S njr = S ra n D n , r . Of course, if S„ ir = 0, set X* r = 0. 
Then, \{D n ^ r ) can be estimated by X* /k n . This leads to an estimate of aj i fc n defined as: 

hkn = Y e i( x r)-jT-, l<i<K, (13) 

r=l Kn 

and consequently to an estimate f n (x) of / n (s) via: 

fn(x) = Y ai,kn e ii x ) = 51 K n (x r ,x)—^. (14) 

i<hn r=l Kn 

Two remarks can be made. First, the estimator does not require knowledge of c, which ensures 
a wide range of applications. Second, the estimator is written as a linear combination of the 
maxima X* r involving the Dirichlet's kernel. Thus, the analysis of the behavior of f n will rely 
on both studies of the Dirichlet's kernel features and of the maxima's stochastic properties. This 
is the topic of the next section. 

3 Basic results 

3.1 Bounds on the Dirichlet's kernel 

For x G [0, 1] and j G {1, 2, 3} define 

B n ,j(x) = Y \Kn(.x r ,x)\ 3 

\r=l 
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and 

B n ,oo{x) = max \K n (x r , x)\ . 

l<r<k n 

In what follows, some theorems involving conditions on B n j(x) for j £ {1,2,3, 00} are fiven. 
Some of these conditions follow easily from the properties of K n , see [12]. In the following lemma 
two properties which are less straightforward are given for the special case of the trigonometric 
basis. 

Lemma 1 Suppose K n is the Dirichlet's kernel associated to the trigonometric basis, 
(i) Ifh n = o(k n ) then sup B njl (x) = O (k n lnh n ), 

X 

(ii) Ifh n lnh n = o(k n ), then for all x £ [0, 1], B n ^{x) ~ (/c n /i n ) 1/2 . 
The proof is postponed to the Appendix. 



3.2 Maxima stochastic properties 

In the sequel, we write: 

K D n,r) = K,r, min f(x) = m n , r , max f(x) = M n>r . 

Recall that X* /r . is the maximum of the second coordinates of the set of points £ n ,r-- 
Noticing that, for < x < m n ^ r , 

P(Xl r <x) = P(N* n (D n , r \ (J n , r x [0, x})) = 0), (15) 

we easily obtain the distribution function F n , r (x) = P(X*^ r < x) on [0,m rajr ]: 

F n ,r( x ) = ex P 



yx K n A n ^ r j 



(16) 



Straight forward calculations lead to the following expansions for the mathematical expectation 
and the variance of X* r , where the knowledge of the C 1 -regularity of / compensates for the lack 
of a precise expression for F n ^ r on ]m n . r ,M n , r [. 

Lemma 2 Suppose f is a function of class C 1 , n = o (k^) and k n = o (n/lnn). Then, 

(i) E(X* ) = k n \ n>r - — + o ^ " 

nc 



i,3 ' 



(ii) Var(Xl 



k 2 

r ) 9 9 



The proof of the following lemma, which is more difficult, is postponed to the appendix. 

Lemma 3 Suppose f is a function of class C 1 and k n = o (n/lnn). Let (t n>r ) be a sequence such 
that t n:r = o(n/k n ) and t n>r = o(k^/n). Then, the characteristic function of (X* — m n ^ r ) can 
be written at point t n>r as: 

1 + itn,T^ c F n A m n,r) + O (|t„, r | ^) + O (n~ s ) 
<t>n,r{tn,r) = 1 , ~ jjT 



1 _i_ 7'+ ^ 



<■>' nc 



with s arbitrary large, and F n ^ r = 1 — F n>r . 
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4 Estimate convergences 

We refer to [12] for a careful study of the bias convergences || /„ — / || 2 — > and || f n — f — > 0. 
Then, we just have to consider (f n — f n ). A complete investigation for mean integrated conver- 
gence, mean uniform convergence, L 2 -almost complete convergence and uniform almost complete 
convergence is available at our website. In order to avoid lengthy developments, we just propose 
the following basic result. 

Theorem 1 Suppose f is a C 1 function. If || f n — / || 2 = o(l) and if sup B rh i(x) = o (n 2 /k n ^j, 
then E(\\ /„-/|| 2 ) = o(l). 

Proof : Introducing the random variable Y n , r = (X* r /k n ) — \ n ,r-> we have 

E(|| /„-/„ || 2 ) = E^Y n , r K n (x r ,Xs)Y n ^, (17) 

< 2^|K n (x r ,x,)|(E(y n 2 r ) + E(y n 2 J) (18) 

r,s 

< 4supi? n , 1 (x)^E(y n 2 r ), (19) 

X , 

r=l 

and the result follows from Lemma 2. ■ 

Corollary 1 If K n is the Dirichlet's kernel associated to the trigonometric basis and f is a 
C 1 function then h n In 1 / 2 h n = o{k n ) and k n (\nh n ) l l 2 = o(n) are sufficient conditions for 
E(\\f n -f\\ 2 ) = o(l). 

Proof : From [12], Proposition 3, /i ra (ln /i n ) 1//2 = o(k n ) entails || f n — / || 2 = o(l). Moreover, 
from Lemma l(i), if k n (ln hn) 1 / 2 = o(n), then 

sup B nA (x) = O (kn In hn) = o (n 2 /k n ) , (20) 
and the conclusion follows. ■ 



5 Asymptotic distribution 

In this section, we present a limit theorem for the distribution of (/„ — E/ n ). A similar result is 
not available for (/„ — /) without reducing the bias, which is done in the next section. 

Theorem 2 Suppose that f is a function of class C 1 . If k n = o(n/lnn), n = o (kn" 2 ^ and 

B n ,oo(x) = o(B n ^{x)), then (nc/ B n ^{x))(f n (x) — E(f n (x))) converges in distribution to a stan- 
dard Gaussian variable for all x G [0, 1] . 



Proof : Denote a n . 



m. 



k n j (nc) and for all x £ [0, 1] introduce ip n:X the characteristic 



function of (nc/ B n ^(x))(f n (x) — E(/ n (x))). It expands as: 



1pn,x( t ) 



exp it 



nc 



k n B n> 2(x) 



5^X n (x r ,x)(a n , r -E(X* r )) 



r=l 



x E 



exp it 



nc 



k n B n)2 (x) 



J2 K n(Xr,x)(X ntr - a n ,r) 



r=l 



(21) 
(22) 
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Consider first the argument of (21): 



T n (x) 



nc 



Kn{x r , x){a n , r - E(X* r )), 



k n B ny2 {x) ^ 

and show it converges to 0. By the Cauchy-Schwartz inequality: 



, m , . , nc B n 1 (x) 



a n>r - E(X* r ) 



nc 

< —777 max 



k 



1/2 



Now Lemma 2 entails 



(23) 



(24) 



n 



r n (x)| < —r-jn ( max(A; n A„ ir - m v ) + o 



,3/2 I ' 



which converges to 0. Thus, introducing i„ jr = n ^ K g^(^ t , we obtain: 



VVi,x(*) ~ expzt ( - — J2 K n( x r,x) ] E 



exp it 



nc 



k n B n ,2{x) fr[ 



Y^K n (x r ,x)(Xl r 



(25) 



= exp it — — ^ K n (x r ,x) Yl </>ra,r(*r, 



B n ,2(x) ^ 



(26) 



r=l 



To apply Lemma 3, we have to verify that t n ^ r = o (n/k n ) and t n ^ r = o (A^/n). The first condition 
is satisfied since 

km i 

nc 



Bn,oo \X) 
B n ,2( X ) 



= o(l). 



The second condition is satisfied as well: 



n 



tn,r " 



n 



n 
kt 



= o(l). 



The characteristic function can be seen to be of the order: 



exp it 



1 



B n ,2{x) ^ 



r=l 



fcn / A; - / n \ 



Consider first the logarithm of term (29). A second-order Taylor expansion yields 



4 1} (*) = £ 



r=l 



it, 



k. 



nc 



In 1 + it n , 



k. 



nc 



O 



' Bj 3 (xy 

B Ux). 



(27) 



(28) 



(29) 



(30) 



(31) 



Since B% 3 (x)/B% 2 (x) < B nj00 (x) / B n>2 {x) = o(l), it follows that Jn\x) ->■ -t 2 /2 as n ->• oo. 
Finally, consider the logarithm of (30): 



> / n 2) = £ ln (! + V)i with = ttn,r^Fn,r{ m n,r) + O ( |i n ,r| ^ ) + O (n S ) . (32) 



r=l 



nc 
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Observe that max r \u n r \ converges to from (27) and (28). Thus, for n large enough \u n ^\ < 1/2 
uniformly in r and the classical identity |ln(l + u n>r )\ < \u n>r \ yields 



, B n i (x) nc ( B n i(x) \ n 2 ,, , . 

B n ,2{x)kl \B ny2 {x) k n 



j(2) 

Therefore, the Cauchy-Schwartz inequality leads to J n 2 ^ — >■ and ip njX (t) — >■ e~* 2 / 2 as n — >■ oo. ■ 

Corollary 2 Suppose K n is the Dirichlet's kernel associated to the trigonometric basis and f 
is a C 1 function. If h n = o(k n ), k n = o(n/lnn) and n = o then, for all x £ [0,1], 

nc{h n k n )~ l l 2 {f n (x) — E{f n {x))) converges in distribution to a standard Gaussian variable. 
Proof : From (7), B n,oo{%) < || K n 1 1 00 — 1 + h n , and from Lemma l(ii), B n p{x) ~ (h n k n ) ' so 

that B ni00 (x)/B ni2 (x) = O ((/l n /A:n) 1/2 ) = O (1) . ■ 

Possible choices of fc„ and h n sequences in Corollary 2 are k n = n 2//3 (lnn) e and h n = (lnn) e for 
e > arbitrary small. These choices entail n(h n k n )~ 1 / 2 = n 2 / 3 (lnn)~ e . 



6 Reducing the bias 

The bias can be decomposed as follows: 

E(/ n - /) = (E/ n - f n ) + (/„ - /), (34) 

where the first term in the sum is the statistical part of the bias, and the second term is the 
systematic part of the bias. First, consider the irreductible bias (f n — /). In order to obtain a 
limit distribution for {nc/ B n ^{x)){f n {x) — f(x)), we need to satisfy the condition 

fe^ (A(I »- /(l) »-°- (35) 

Introduce S n (f) = ^Ojej. Equation (3.12) in [12] provides sharp bounds for (S n (f) — f n ), so 
i=o 

that the question reduces to considering (S n (f) — /), which only depends on the basis. We shall 
see that the trigonometric basis satisfies (35) under reasonable conditions on h n and k n . In the 
case of a general C 1 basis, we shall take (35) as a condition. 
Now, it follows from Lemma 2(i) that 

E fn{x) ~ fn{x) = ]T K n ( Xr ,x) (^Ihll - \\ ( 36 ) 

presents a negative component which should be eliminated. To this end, for r = 1, . . . , k n , define 
Z* r by Z* r = if S Ttjr = and Z* r is the infimum of the second coordinates of the points of 
Y> n , r otherwise. Then, the random variable 

Z n = 7T X! Z n,ri (37) 



has a mathematical expectation 



E(Z„) = ^+of^), (:)8, 



nc \k n 
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and a variance 



Var(Z n ) 



n 2 c 2 



(use Lemma 2 for Z* r ). We define a corrected estimate by 

fn{x) = jrK n (x r ,x) ( "' r , — - ) = /„(x) +ff n (x). 



r=l 



Lemma 4 Suppose f is a function of class C 1 , n = o (k 2 ) and k n = o (n/lnn). T/ien, 
— ^— £(/ n (x)) - /n(x) = o (n 2 /k 7 J 2 ) , Vx G [0, 1]. 

//, moreover, for all x G [0, 1], A; n = o (-B 2 2 (x)^ and /35j holds, then 



nc 



Proof : We have 
nc 



B n ,2(x) 



B n;2 (x) 

E(/n(x)) - f n (x 



(g n (x) - E{g n {x))) = o P (l). 



nc.B n i(x) 
< — max 

«n %,2 I l<r<fe„ 



E(X* r ) + E(Z n ) - fc„A n ,, 



nc I!,,., (x) / n \ / /r 



&n B n>2 {x)° \k n 



1,7/4 I ' 

, "'71 / 



from (38), Lemma 2 and the Cauchy-Schwartz inequality. 
Now, applying (35) to the constant function / = 1 yields 

— ^2K n (x r ,x) 1, 



r=l 



as n — > oo. Therefore, Var(<7„(x)) ~ Var(Z n ) ~ k n /(n c ) and then 



Var 



nc 



#n,2(x) 



(<?n(x) - EQ n (x)) 



kn 



Bl, 2 {x) 



which converges to 0. 



(39) 



(40) 



(41) 
(42) 



(43) 



(44) 



Theorem 3 Suppose that f is a function of class C 1 . If the following conditions are verified 
(i) k n = o (n/lnn), n = o (kf/ 2 ^j, 

(ii) for all x G [0,1], max(/c n /2 , B n>00 (x)) = o(B n>2 (x)), 



(Hi) for all x G [0, 1], 



nc 



B n ,2( x ) 



\f n (x)-f(x)\=o(l), 



then, for allx G [0, 1], (nc/ B nt2 (x))(f n (x)— f (x)) converges in distribution to a standard Gaussian 
variable. 

The proof is a simple consequence of the expansion 

(fn ~f) = (fn ~ E(/ n )) + (g n - E(g n )) + (E(f n ) - f n ) + (f n - /), (45) 
and of Theorem 2 and Lemma 4. 
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Corollary 3 Suppose K n is the Dirichlet's kernel associated to the trigonometric basis and f is 
a C 2 function. If h n lnh n = o (k n ), n = o (hf/ 2 kj/ 2 ^J , nhl/ 2 In h n = o (kn 2 ^j and k n = o (n/ln n) 

then, for all x € [0, 1], nc(h n k n )~ l l 2 (f n (x) — f{x)) converges in distribution to a standard Gaus- 
sian variable. 

Proof : Conditions (i), (ii) of Theorem 3 are verified. Consider (Hi). On using (8) we have 

nc(h n k n )~ 1/2 \S n (f)(x) - f(x)\ < nc(h n k n )~ 1/2 ^2 \ai\\ei(x)\ 

i>h n 

< V2nc(h n k n )~ 1/2 ^2 \ a i\ 



i>h n 



< V2nc(h n k n )- 1/2 J2 *" 

i>h n 



A straightforward calculation yields 

nc 



B n>2 {x) 

From [12], equations (3.11) and (3.12) 



S n (f)(x)-f(x)\ = 0(nh- 3 / 2 k- 1 



/2 



\ s n(f)(y)-fn(y)\ = o^j o 1 



dK, 



dx 



-(v,y) 



dv 



and 



Therefore, 



/' 

Jo 



dK n 



dx 



dv = O (h n In h n ) 



nc 



\S n (f)(y) - fn(y)\ = O (nk-V 2 h}/ 2 lnh r 



(46) 
(47) 

(48) 

(49) 
(50) 



and (47) with (50) conclude the proof. ■ 

Possible choices of k n and h n sequences in Corollary 3 are k n = n 4//5 (lnn) 3//5 (lnlnn) e and 

h n = n 2//5 (lnn) _1 / 5 (lnlnn) e for e > arbitrary small. These choices entail n(h n k n )~ l l 2 = 
n 2//5 (In n) _1 / 5 (In In n) _£ . 



7 Conclusion and further developments 

In this paper, we showed how the convergence results established in [7] in the case of the Haar 
basis can be adapted for any C 1 basis under some assumptions on the Dirichlet's kernel behavior. 
We have emphasized that the estimator and these assumptions only depend on the Dirichlet's 
kernel of the basis. This suggests to define a new estimator based on a Parzen-Rosenblatt kernel. 
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Appendix 

Proof of Lemma 1 

(i) In the sequel, we shall use the inequality found in [12], equation (2.11): 

s'm(pu) 



sin(u) 



7T 



< pl[ 0l i](|«|) + ^-rrl[5,7r/2](kl), 



(51) 



for all p > 0, < 5 < tt/2 and \u\ < tt/2. Taking account of the periodicity and symmetry 
properties of the trigonometric kernel, it suffices to study 

2 [W2] + l 

sup — \K n (x r ,x)\, (52) 

x6[0,l/fc„] K n r=1 



where [u] denotes the integer part of u. Let us write 

1 [k n /2]+l 1 [ 7 „] 1 [fc„/2]+l 

^ ' K n (x r , x) = -j— ^ ' K n (x r , x) + ^ ] K n (x r , x), 



A: 



(53) 



r=l r=l ■"» r=[ 7 n]+l 

with 7„ = k n /(h n + 1), and consider the two terms separately. 



Introduce 5 = (^/k n )([y n ] — 1/2). For r = 1, . . . , [y n ], we have 7r(x r — x) < 5 and thus 
(51) yields \K n (x r , x)\ <l + h n which gives in turn 



[in 



^2 \K n (x r ,x)\ < 1. 



(54) 



r=l 



• For r = [j n ] + 1, . . . , [k n /2] + 1, we have ir(x r — x) > 5 and consequently (51) yields 

[fen/2] 



1 [fe„/2]+l [fe„/2]+l 

— 53 |K„(x r ,x)i<— 53 



r=[7n]+l 

Therefore, 



2fc„ ^ l 



fe « r = [7n]+ i 2 ^ x ) 2 fc ™ ,= [7n] £ (r - i) 



[fe„/2]+l + J i 

f £ |Jr„(* r ,*)|<i/; 2fc "^< ±in(4(*» + i)), 



(55) 



(56) 



for fc n > 2(/i„ + 1). 
Finally, collecting (54) and (56), we obtain 

1 



SU P TrJ2\ K n( x r,x)\ < 2 + ln(4(/i n + l)). 



xe[0,l] r=1 



(ii) From [12], equation (4.14), 

? 2 

5 n,2 



£^(x) 



- 1 



< II Kn Hoo 



K n (x, x) fe 



k n K n (x, x) 

In the case of the trigonometric basis K n (x, x) = \\ 
equation (3.12) we have 



1 f 1 dK 
-sup / -——(x,v) 

In x Jo dy 



dv. 



(57) 



(58) 



l + h n (see (7)) and from [12], 



Jo 



dy 



■(x,v) 



dv = (h n In h n ) = o (k n ) . 



(59) 



The conclusion follows from (58) and (59). 
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Proof of Lemma 3 

Consider the expansion 



E(e 



P(X* r = 0) + 



Jn 



M n , r 



e lxt ^F n)r {dx). (60) 



The first and second term can be computed explicitely since (16) provides an expression of F n ^ r 
on [0,m njr ]: 



E(e 



itn,rX*, : 



+ ! ^n,rtm TO , r j e / e ixt n ,r Fnr r dx) . (61) 

^ Tfln.r 



1 + itn,r nc 



The third term can be expanded as 



/ e^'-F^dx) = e^^F^m^) + / (e" 4 "' - e lt ^ m ^)F n ^ T {d. 



with |e l:r * n ' r - e**".'-" 1 ".'- 1 < (M n>r - m„ ir ) |i n ,,-|- Thus 

pM„, r 



(e 



)F n , r (dx) 



< {M n ^ r tnn,r) \tn,r\ F n ^ r (jn n ^ r ) 



nc 



^ (-^n,r m n,r) \tn,r\ , (^n^ 



n,r m n,r 



n 



- [KA k 3.) ■ 

Collecting (61)-(63), and remarking that k n = o(n/\nn) yields 
the proof. 



(62) 



(63) 

= O (n~ s ) concludes 
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